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TECHNICAL NOTE 2826 

SIMULATION OF LINEARIZED DYNAMICS OF GAS-TURBINE ENGINES 
By J. R. Ketchum and R. T. Craig 


SUMMARY 

Through the use of an electronic analog computer in the simulation 
of controlled aircraft -engine performance, one method of engine simu- 
lation has proved to he preferable from many considerations . 

The equations used in developing this method of simulating the 
dynamics of gas-turbine engines are derived in general form from engine 
functional relations. This general simulation method can be utilized 
in the consideration of any first-order linear system and is designed 
for use in conjunction with control components for small perturbation 
or stability studies of controlled system operation. 

A simulation of the response of a turbojet engine to a step change 
in an independent variable is made, and comparison of the experimental 
and simulated results indicates the validity of the simulation method 
presented. 

Limitations on the use of altitude and flight-speed generalization 
factors in determining the coefficients necessary for the simulation of 
engine dynamics axe discussed. 


INTRODUCTION 

A considerable amount of present-day control- system design and 
analysis is accomplished through the use of simulation techniques . 

With such techniques, part or all of a physical system is replaced by 
its mathematical representation, usually in the form of an analog com- 
puter facility. The mathematical representation then is used in sub- 
sequent examinations of controlled system behavior (reference l) . 

The first step in the process of simulation is the determination 
of equations descriptive of the behavior of the system under all con- 
ditions of operation. The second step requires that these equations 
be put into a form applicable to the computer or simulation facilities 
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available. Although the details of setting up a problem will vary con- 
siderably because of the wide variety of such machines, the same general 
equations are used in all cases. 

An electronic analog computer has been used for some time at the 
NACA Lewis laboratory for simulation of aircraft -engine performance and 
studies of controlled engine operation. During the course of these 
activities, one method of simulation has proved to be the most advan- 
tageous in economy of computer facilities, utilization of experimental 
data, and accuracy of results. The object of this paper is to develop 
and present this method and to show its applicability to several types 
of gas-turbine engine. 


BASIC DYNAMICS OF GAS-TURBINE ENGINES 

For some given equilibrium condition, there will be established 
definite values of engine rotational speed, temperatures, pressures, 
and torques. At an equilibrium condition, the torque produced 
by the turbine will be totally absorbed by the compressor. If 
the equilibrium is disturbed by the injection of additional fuel, 
the instantaneous torque produced by the turbine will be 
greater than that absorbed by the compressor, and the engine will 
accelerate to a new equilibrium condition. For most practical cases, 
the rate of speed increase is proportional to the difference in tur- 
bine and compressor torques. This condition is typical of first-order 
systems. Consequently, exponential equations may be employed to des- 
cribe the dynamic operation. 

For description of the dynamic operation, a term is required which 
will determine the time history of the engine responses. In view of 
the exponential nature of the engine behavior, these responses can be 
described by specifying a time constant. This time constant is a 
ratio between the energy storage and energy dissipative elements and 
is the time required for 63 percent of a change to occur. It is a 
transient term and requires transient data for its evaluation. 

The extent to which a variation in fuel flow will change 
an output variable such as engine rotational speed is the 
gain, or sensitivity of the engine rotational speed to changes in 
fuel flow. At any operating point, this gain is evaluated as the slope 
of the engine speed - fuel-flow curve and is the partial derivative 
dN„ 

— — evaluated with all independent variables other than fuel flow held 

dW 

constant. In a particular engine there are gain terms relating changes 
of rotational speed to changes in fuel flow, exhaust -nozzle area, and 
propeller blade angle; there axe also gains which relate changes in 
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pressures, temperatures, and. "torques "to "these same variations in fuel 
flow, exhaust-nozzle area, and blade angle. Gains related to changes 
in inlet geometry, guide-vane angle, and other independent variables 
also may be considered. The number of these terms which must be con- 
sidered will depend upon the type of engine and the particular analy- 
sis or study undertaken. The gain terms will determine the quantita- 
tive change in a dependent or output variable produced by the variation 
of an independent or input variable. It will give no idea of the rate 
at which this change occurs. 

To obtain the time relation of the various parameters, certain 
partial derivatives, evaluated with all independent variables held 
constant, must be obtained. These terms generally require the use of 
transient data for their determination, and they are specifically 
referred to in later sections of this report. All the gain terms and 
the time constant can be treated with standard altitude and flight - 
speed generalization factors . 

Linear analysis was used in the development of this method of 
simulation. Although some broad simplifying assumptions must be made 
to allow use of linear analysis, it is possible to obtain reasonably 
accurate results with a minimum of mathematical complications . It has 
been demonstrated by comparison of analytic results with experimental 
performance that linear analysis is adequate to permit prediction of 
stability of engine and controls systems and to predict the response 
of these systems to small perturbations (reference 2). An example of 
particular importance in the use of linear analysis is the study of 
afterburner effects on the primary controlled engine. 


DESCRIPTION OF ANALOG COMPUTER 

The method of simulation of gas-turbine engines that is discussed 
in this report was developed for use on the electronic analog computer 
at the NACA Lewis laboratory. This computer is composed of operational 
amplifiers as indicated in figure 1. With a transfer function defined 
as the ratio of the Laplace transform of the output variable to the 
Laplace transform of the input variable, the approximate transfer func- 
tions of the operational amplifier is (reference 3) 

y(s) = _ Z 2( s ) 
x(s) Z 1 (s) 

(Symbols are defined in appendix A.) By selection of suitable pure 
resistive impedances, amplifiers may be constructed which will change 
algebraic sign, multiply by constant -valued coefficients, add, and sub- 
tract, If a capacitor Cg is utilized for impedance Zg and a 
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resistive element Rq for impedance , the amplifier unit yields as 
an output a time integral of the input variable. The resulting unit 
will have a time constant equal to RqCg. A number of computational 
elements have been assembled to form the complete computer shown in 
figure 2 . 

Of special interest in connection with this report is the so-called 
matrix unit, composed of sign- changing, coefficient, and summation units, 
which yields the following relation between input and output variables: 


y x = a 11 x 1 + a 12 x 2 + a 13 x 3 + a 14 x 4 

^2 = a 21 x l + a 22 x 2 + a 23 x 3 + a 24 x 4 
^3 = a 31 x l + a 32 x 2 + a 33 x 3 + a 34 x 4 
y 4 = a 41 x l + a 42 x 2 + a 43 x 3 + a 44 x 4 


A diagrammatic representation of this unit is shown in figure 3. 

Time behavior of variables is displayed on oscilloscopes and can 
be recorded by photographic techniques. 


DEVELOPMENT OF SIMULATION METHOD 

Turbine -propeller -engine simulation . - Development of the equa- 
tions describing dynamic response of gas-turbine engines is based upon 
the assumption that thermodynamic and flow processes are quasi-static; 
that is, these processes continuously progress from one state of equi- 
librium to another on an equilibrium curve (reference 4) . This assump- 
tion permits the writing of functional relations between the various 
input and output variables . 

In analysis of turbine -propeller engines, the most useful rela- 
tions are those between engine rotational speed, turbine-outlet tem- 
perature, fuel flow, and propeller blade angle. For analytic purposes, 
the accelerating torque is considered also. This torque is the differ- 
ence between the torque produced by the engine and that absorbed by the 
propeller. The engine torque is a function of fuel flow and engine 
rotational speed: 


% = q(w,» e ) (i) 

The propeller torque is a function of propeller blade angle and propel- 
ler rotational speed: 
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<Jp = f 2 (^V (2) 

The unbalanced torque referred to the engine shaft is the difference 

between engine and propeller torques, (Q e - — ) , where r is the gear 
ratio . 

The turbine -out let temperature is taken to be a function of fuel 
flow and engine rotational speed: 

T = f 3 (W,N e ) (3) 


In a similar manner, any other engine-output variable may be 
expressed as a function of fuel flow and engine rotational speed. 


Engine rotational speed is related to the unbalanced torque by the 
expression 



(4) 


where 1^ is the sum of the polar moment of inertia of the engine and 
the polar moment of inertia of the propeller referred to the engine. 

Functions (l) and (2) may be expanded and linearized to give 




dW 


AW 


N c 


+ 


dQe 

SN e 


AN e 

W 


(5) 



dQp 

S|T 




AN p 

P 


( 6 ) 


where 

dW 

to fuel flow, with engine speed held constant . Other terms axe to be 
interpreted in similar manner. 

The accelerating torque is then 


is the partial derivative of engine torque with respect 
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AQ = AQ e - AQ p 


SQe 

Sw 


l 

AW - - — ^ 


N. 


r Sp 


AP + 


JL 
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r OH p 


AN 


( 7 ) 


Equation (3) is expanded and linearized to yield 


AT = 


dT 

dW 


AW + 


N, 


St 

dKL 


AN C 


w 


( 8 ) 


The values of the terms forming the coefficients may he found 
from the slopes of steady-state operating curves, evaluated at par- 
ticular values of fuel flow, blade angle, engine rotational speed, 
altitude, and flight speed. These equations together with equation (4) 
are introduced into the computer using the matrix unit described in 
a previous section. The interconnection diagram is shown in fig- 
ure 4(a). This simulation will give an adequate description of dynamic 
operation of turbine-propeller engines over considerable ranges of 
operation (reference 5) . 

Turbo jet -engine simulation . - Development of equations describing 
dynamic operation of turbojet engines requires the basic assumption of 
quasi-static operation, just as required in the case of turbine- 
propeller engines . 

The input variables of turbojet engines are fuel flow and exhaust- 
nozzle area. The output variables generally considered in controls 
analysis are engine rotational speed and turbine -out let temperature. 
Another variable is necessary to account for engine operation at non- 
equilibrium conditions . The variable considered is accelerating 
torque, which is the difference between the torque produced by the 
turbine and that absorbed by the compressor. 

The relations assumed between the engine variables are 

Q = f 3 (W,A,N e ) (9) 


and 


T = f 4 (W,A,Q) 


( 10 ) 
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Under steady-state conditions, unbalanced torque is zero and engine 
rotational speed is an implicit function of fuel flow and exhaust- 
nozzle area. These functions may be expanded and linearized about 
some equilibrium point to yield 


d W 


A,N, 


AW + |» 

dA 


AA + 


W,N P 


5il 


AIL 


W,A 


5t 

AW 

+ 

AA 

+ 

5w 

A,Q 

3a 

W,Q 

5q 


AQ 


W, A 


( 11 ) 


( 12 ) 


To evaluate these coefficients in this form would require plots of 
accelerating torque as a function of fuel flow, engine rotational speed 
and temperature as functions of accelerating torque, and temperature 
as a function of fuel flow. There is, however, no direct way of meas- 
uring this torque in a turbojet engine; consequently, to obtain the 
required plot would necessitate either a considerable amount of data 
processing or impractical test procedures . It is possible to circum- 
vent these difficulties by manipulating the torque equation (ll) in 
the following manner: 


But 


and 


AQ 



dQ 

5w 


A,N e 


5q 


AW 


W,A 


5Q 


5a 

W,N e 

5q 


5n 0 

W,A 
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5Q 
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Therefore , 


AQ 


Sq 

Sn, 


W, A 



+ 


c®e 

SA 


AA - AN e 
W 


(13) 


This manipulation has yielded terms whose coefficients may he evaluated 
from the slopes of the steady-state plots of speed against fuel flow 
and area. The slopes are evaluated at the appropriate operating points. 

Sq 

Equation (13) is multiplied by 

1N e 

tion (12) in this manner: 



and substituted in equa- 

W.A 


But 



Furthermore, in steady state, Q - 0, hence 


ST 

Sw 


A,Q 



♦ 


and 


ST 

SA 


W,Q 


St 

Sa 
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Both of these terms are obtainable from steady-state data. 

Thus, the equation of the response of turbine -outlet temperature 
is 
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(14) 


The torque-speed relations of a turbojet engine also may be 
expressed by equation (4) . This equation is multiplied and divided by 
SQ 

in the fo Honing manner: 

W,A 





The engine time constant is defined 


9 dQ 



W,A 


in the preceding equation to yield the expression 


and substituted 



(15) 


St 

The engine time constant, as well as — 

oN fi 


, is a typically 


W,A 


transient term and requires transient tests for its evaluation. 


Equations (13) to (15) will determine speed and temperature 
response to changes in fuel flow. These are the equations used in 
simulation of turbojet engines. An interconnection diagram showing 
a method of using computer components for solution of these equations 
is shown in figure 4(b). (Terms in any column operate on the respective 
input to that column, and outputs are summed across the rows.) 

Simulation of additional variables . - If additional input varia- 
bles are to be considered, equations (9) and (lO) must be extended in 
this manner: 
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Q - f^(W, A,N e , . . • X) 

T = f g (W,A,Q, . . . X) 

where X denotes any additional input variable, such as variable inlet 
diffuser geometry, for example. 

Functional relations of additional output variables can be formed 
by the same basic reasoning as that used for the temperature relations . 

If any additional output variable (compressor-discharge pressure, for 
example) is denoted as G, then 

G = f 2 (W,A,Q, . . . X) 

which, upon suitable expansion, becomes similar in form to equation (14) . 

The manner in which the basic method may be extended to include 
simulation of additional engine variables in orderly fashion is illus- 
trated in figure 5 (reference 6) . 

The manner in which the equations are evolved results in an array 
of equations wherein a maximum number of terms can be evaluated from 
steady-state data. All coefficients are slopes of steady-state curves, 
except the engine time constant x e and the terms in the right column 
(fig. 5) . Methods of determining the required values of terms with the 
minimum amount of testing are discussed in reference 7 . 


SAMPLE SIMULATION 

The method shown in figure 4(b) was used to simulate the response 
of turbine -outlet temperature and engine rotational speed of a turbojet 
engine to an approximate step change in fuel flow. The simulated 
responses are shown in figure 6, and the comparable experimental responses 
are shown in figure 7. In this engine-dynamics simulation it was necessary 
to include dynamics of fuel system and test facilities . The fuel sys- 
tem was simulated by two first-order time lags and an interposed clip- 
ping network whose response duplicated the time behavior of fuel flow; 
the speed and temperature sensors were simulated by first-order time 
lags . The altitude test facilities were such that during transient 
operation engine-inlet pressure variations could not be eliminated, 
and thus a variation in nominal flight speed was produced. Compensa- 
tion for this extraneous effect was introduced by evaluating gains and 
time constants directly from oscillograms that were to be simulated. 

In figure 8 the simulated and experimental responses are compared. 

The degree of correlation indicates the validity of the description 
and simulation of the engine dynamic behavior. 
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All experimental data obtained to date indicate that turbojet 
engines remain first order and linear for small perturbations over a 
wide range of flight conditions. Therefore, the developed method of 
simulation permits the exploration of a wide region of engine opera- 
tion, provided only that the coefficients of the equations can be 
evaluated. If a complete set of equation coefficients can be obtained 
for any one flight condition, the coefficients may be adjusted for 
use at other flight conditions through the use of standard generaliza- 
tion factors. Complete generalization of engine dynamics should not 
be done indiscriminately, however. The fact that some knowledge must 
first be obtained of the altitude effects on the particular engine 
considered is discussed in appendix B. 


CONCLUDING REMARKS 

An electronic analog computer has been used for some time at the 
NACA Lewis laboratory in simulation of controlled aircraft -engine per- 
formance . One method of engine simulation has emerged as the most 
advantageous in economizing computer facilities, requiring the mini- 
mum amount of experimental data, and attaining the most accurate 
results . 

The equations used in developing this method of simulating the 
dynamics of gas-turbine engines are derived in general form from engine 
functional relations . This general simulation method can be utilized 
in the consideration of any first-order system that can be linearized 
over a portion of its operating range. Any number of dependent and 
independent variables can be included in the analysis and simulation. 

A simulation of the response of a turbojet engine to an approxi- 
mate step change in fuel flow is made, and comparison of the simulated 
and experimental results indicates the validity of the simulation 
method. Since gas-turbine engines generally remain first order and 
linear for small perturbations over large operational regions, the 
simulation method can be utilized to explore engine characteristics 
over any range of flight conditions . 

The use of altitude and flight-speed generalization factors in 
determining the equation coefficients necessary for the simulation of 
engine dynamics is discussed, and attention is drawn to the limita- 
tions placed on this use by the nature of component efficiency varia- 
tion with flight condition. 


Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio 
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APPENDIX A 
SYMBOLS 

The following symbols are used in this report: 

A exhaust-nozzle area 

a constant (coefficient) 

C capacitance 

G additional output variable 

I polar moment of inertia 

K amplifier gain 

N rotational speed 

Q unbalanced (accelerating) torque 

R resistance 

r gear ratio, engine speed/propeller speed 
s complex Laplacian operator 

T absolute turbine-discharge temperature 

t time 

W fuel flow 

X additional input variable 

x arbitrary input signal 

y arbitrary output signal 

Z impedance 

(3 propeller blade angle 

A deviation from initial operating point 

combustion efficiency 
1 time constant 
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Subscripts : 

• 

e 

engine 


f 

final 


i 

initial 

2683 

P 

propeller 

t 

total 
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APPENDIX B 

USE OF GENERALIZATION FACTORS IN EVALUATION 
OF ENGINE DYNAMIC CHARACTERISTICS 

The developed method of simulation permits the exploration of engine 
operation over a wide range of flight conditions, provided only that the 
coefficients of the equations can he evaluated. 

Standard generalization factors have found widespread use in the 
handling of engine steady-state operational characteristics, since they 
allow prediction of engine performance in regions not covered by test. 

A similar use can be made of the generalization factors in the handling 
of dynamic characteristics. In the discussion of generalization factors 
it is noted that although many performance maps will generalize, those 
directly involving component efficiencies which vary greatly with flight 
condition will not do so . Fuel-flow relations are of the type having 
these large variations, and fuel is of primary interest with respect to 
engine controls. The generalization of engine dynamic terms, however, 
is not as limited as is the generalization of static characteristics . 

In some cases it has been illustrated that even though the absolute val- 
ues of the fuel-flow relations would not generalize, the slopes of the rela- 
tions (the dynamic characteristics) would do so. Such action indicates that 
even though the values or levels of efficiencies varied as flight conditions 
were changed, the form of the variation of efficiencies with rotational 
speed remained constant . 

Some limitations must be placed on this use of generalization fac- 
tors, however. For the engine with the dynamic characteristics shown 
in figure 7, it was found that the gain terms such as Sn/^W would not 
generalize. Figure 9 illustrates that the generalized gain terms varied 
as the altitude of operation changed. This variation indicates that not 
only efficiency levels but also the effect of rotational speed on com- 
ponent efficiencies changed with altitude. 

The variation of combustion efficiency with flight condition, engine 
speed, and exhaust-nozzle area was determined} and compensation for this 
variation was added to the data previously presented in figure 9. Fig- 
ure 10 shows that the gain terms of speed - fuel flow now will 
generalize. Other generalized gain terms (dT/dW, for example) show 
similar improvement when compensation for combustion efficiency varia- 
tions is introduced. 

For the particular engine under consideration it appeared that 
efficiency changes in other components such as the compressor and tur- 
bine did not have a noticeable effect on the engine gain terms. The 
uncompensated gains could not be considered to generalize, but compensa- 
tion for combustion efficiency variations alone was sufficient to obtain 
accurate results in determining the necessary gain terms of a turbojet 
engine at altitude. 
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It therefore may he concluded that complete generalization of engine 
dynamics should not he done indiscriminately. Some knowledge must first 
he obtained of the altitude effects on the particular engine considered. 
If efficiencies vary consistently with altitude as did those considered 
in reference 4, complete generalization of dynamic terms may he used. 

For other engines sufficient altitude test data must first he obtained 
and the variations of efficiencies incorporated into the dynamic terms . 
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Figure 2. - Electronic analog computer used for simulation of dynamics of gas-turbine engines and investigation of 

controlled engine operation. 
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Figure 3. - Schematic diagram of matrix unit of electronic 

analog computer. 
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(a) Turbine-propeller engine. 
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(b) Turbojet engine. 

Figure 4. - Interconnection diagram used in simulation of linearized dynamics 
of jet engines with electronic analog computer. 
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Figure 5. - Interconnection diagram used in simulation of linearized dynamics of 
turbojet engine showing manner in which additional variables can be simulated. 
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Figure 6. 
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- Simulated speed and temperature responses of turbojet engine to change in fuel 

flow. 
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Figure 7. - Experimental response of turbojet engine to change in fuel flow at constant area and flight speed. 
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Figure 8. - Comparison of simulated and experimental responses of turbojet engine. 



Figure 9. - Variation of generalized gain with engine rotational speed for turbojet engine at 

constant flight speed. 
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Figure 10. - Variation of generalized gain with engine rotational speed for turbojet engine. Gain compensated 

for combustion efficiency variations. 
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